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The warm inflation paradigm considers the continuous production of radiation during inflation 
due to dissipative effects. In its strong dissipation limit, warm inflation gives way to a radiation 
dominated Universe. High scale inflation then yields a high reheating temperature, which then poses 
a severe gravitino overproduction problem for the supersymmetric realisations of warm inflation. In 
this paper we show that in certain class of supersymmetric models the dissipative dynamics of the 
inflaton is such that the field can avoid its complete decay after inflation. In some cases, the residual 
energy density stored in the field oscillations may come to dominate over the radiation bath at a 
later epoch. If the inflaton field finally decays much later than the onset of the matter dominated 
phase, the entropy produced in its decay may be sufficient to counteract the excess of gravitinos 
produced during the last stages of warm inflation. 

PACS numbers: 98.80.Cq, 11.30.Pb, 12.60.Jv 



I. INTRODUCTION 

The flatness required of inflationary potentials, to sat- 
isfy density perturbation constraints, most commonly re- 
lies of Supersymmetry to protect against radiative cor- 
rections. Although SUSY is hugely successful in build- 
ing models of inflation, associated with it are also some 
problems. One of these is the gravitino overproduction 
problem. Very simply the problem is, ending inflation at 
too high a temperature can lead to a gravitino abundance 
that is prohibited by nucleosynthesis constraints when we 
consider a massive unstable gravitino (in this paper we 
will mainly discuss the case of the massive unstable grav- 
itino with the mass of 777,3/2 ~ 0(1) TeV). One alternative 
is to end inflation at a lower temperature, but there are 
also advantages to a high temperature exit from inflation, 
notably it can lead to effective leptogenesis 0, Q . 

There are two distinct dynamical realizations of infla- 
tion, cold and warm inflation. Cold inflation is the stan- 
dard scenario in which the inflaton is assumed to be non- 
interacting during inflation, and thus there is no particle 
production during inflation, so the Universe supercools 
|3[ . Only after inflation, interactions of the inflaton with 
other fields is assumed significant, and a reheating phase 
occurs when the vacuum energy used to drive inflation is 
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converted into particles that forms the subsequent radia- 
tion dominated phase. In the alternative warm inflation 
picture, the inflaton interacts with other fields during the 
inflation phase, which leads to particle production con- 
current with inflationary expansion Q (for recent reviews 
please see (E-01)- The presence of a radiation energy 
density is not inconsistent with the General Relativity 
requirements for realizing inflation, which only requires 
that the vacuum energy dominates the energy density in 
the Universe. Thus in the warm inflation picture, the 
presence of radiation during inflation implies this phase 
smoothly ends into a radiation dominated phase without 
a distinctively separate reheating phase. This offers an 
alternative dynamic solution to the graceful exit prob- 
lem of inflation. The presence of radiation then implies 
the fluctuation of the inflaton, which are the primordial 
seeds of density perturbations, are now thermal 
rather than the quantum fluctuations that occur in cold 
inflation. 

The equations of motion for the inflaton field and for 
the radiation density in the presence of a dissipation 
mechanism are 

4> + 3H(l + Q)cf) + V' = 0, (1) 

and 

p r + 4H Pr = T0 2 , (2) 

where Q = T/3H, T is the dissipative coefficient, over- 
dots stand for time derivative, and ' = Throughout 
the paper we use natural units c = h = fee = 1 and 
Newton's gravitational constant is 8ttG = nip 2 , where 
777p = 2.4 x 10 18 GeV is the reduced Planck mass. 
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Since interactions are important during warm infla- 
tion, if the fields interacting with the inflaton are at high 
temperature, then it is difficult to control the thermal 
loop corrections to the effective potential that is needed 
to maintain the very flat potential required for inflation 
(ill Il2l |. However, if the fields interacting with the in- 
flaton are at low temperature, then supersymmetry can 
be used to cancel the quantum radiative corrections, and 
maintain a very flat inflaton potential. In this paper 
we consider a class of supersymmetric models for which 
the dissipation coefficient T has been computed in the 
equilibrium approach [ll| for the low-temperature regime 
[l3l [l4| . The dissipation mechanism is based on a two- 
stage process [Hj]. The inflaton field couples to heavy 
bosonic fields, x an d fermionic fields ip x , which then de- 
cay to light degrees of freedom. These light degrees of 
freedom thermalize to become radiation. The simplest 
supcrpotential containing such an interaction structure 
is 



W = g<PX 2 + hXY 2 , 



(3) 



where X and Y denote superfields, and <f>, x an( i V 
refer to their bosonic components. Such an interaction 
structure is common in many particle physics SUSY mod- 
els during inflation, the field y and its fermionic partner 
y remain massless or very light, whereas the field x arL d 
its fermion partner ip x obtain their masses through their 
couplings to <p, namely = m x = gcj>. The regime of 
interest is when m x , m,^ > T > H, and this defines what 
is usually referred to as the low-temperature regime. For 
this regime the dissipation coefficient, when the super- 
fields X and Y are singlets, is found to be [TH, fl4| 



0.64 g 2 h 4 
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mi 



(4) 



where T is the temperature of the radiation bath, 
p r = C r T A , where C r = 7r 2 <7*/30 and <?* is the number of 
relativistic degrees of freedom. The above dissipative co- 
efficient is calculated under the adiabatic approximation, 
and associated with it arc consistency conditions which 
ensure that the microscopic dynamics is faster than the 
macroscopic motion [Tl| , 



' X (V>x) 



> r ,H. 



(5) 



In supersymmetric theories C r ~ 70. However, the su- 
perfields X and Y may belong to large representations 
of a GUT group. In that case, the dissipation coefficient 
picks up an extra factor J\f = N x N 2 cc&y , where N x is the 
multiplicity of the X supcrficld and A/d ecay is the num- 
ber of decay channels available in X's decay. Typically, 
having enough dissipation during inflation requires large 
multiplicites, Af x ~ A/decay ~ 0(100). In that respect, 
string theory can be a natural place for warm inflation 
due to the presence of large numbers of moduli fields fl6| . 

Following this approach, it has been recently shown 
[T7I li~8j that chaotic and hybrid inflation models may 



support some 50 to 60 e-foldings of warm inflation in the 
strong dissipative regime, but such models can lead to 
an overproduction of gravitinos (l9j . In the context of 
hilltop models, it was found in Ref. (2(| that although 
warm inflation agrees with current observations the tem- 
perature of the radiation at the end of inflation exceeds 
the current bounds on thermal gravitino production for 
massive gravitino with its mass m 3 / 2 = O(TeV). 

In order to dilute the excess of gravitinos thermally 
produced towards the end of warm inflation, it is possi- 
ble to argue that the necessary entropy production owes 
to the decay of the inflaton field, which comes to domi- 
nate the Universe at a later epoch. However, if the ratio 
Q ^s> 1 , the inflaton field decays completely right after in- 
flation and no later entropy production can be attributed 
to it. It is worth emphasizing though, that T, as given 
in Eq. (|4|, is time-dependent. Moreover, it always falls 
faster than the Hubble parameter during the radiation 
dominated epoch that follows after inflation: T oc a~ 3 
whereas H oc a -2 . Therefore, if Q is not too large, it is 
plausible that soon after inflation the system moves into 
the weak dissipation and hence the inflaton field does not 
decay completely until a later epoch. Once in the weak 
dissipation regime, the average scalar density of the field 
decreases as a -3 , hence it may come to dominate over the 
radiation density thus driving a late matter-dominated 
epoch. This epoch is terminated by the perturbative de- 
cay of the inflaton field. The entropy produced by this 
decay may be sufficient to dilute the excess of gravitinos. 
After this, the inflaton decays completely before the Big- 
bang nucleosynthesis epoch and the thermal bath of the 
hot big bang is recovered 1 . 



II. FIELD EVOLUTION 

In order to investigate the field evolution during infla- 
tion and its subsequent oscillatory phase we consider a 
class of hilltop models with the scalar potential given by 



V = v f(4>) , 



(6) 



where Vb is a characteristic density scale and / is a di- 
mensionless function with a maximum at <j) — and a 
minimum at 4> = 4>v, where it vanishes. We parametrise 
the field's expectation value by the dimcnsionlcss quan- 
tity 8, defined by 



4>=4> v {l-5). 



(7) 



In the following we use subscripts "*" , "e" and "o" to de- 
note the time when cosmological scales exit the horizon, 



1 When we do not consider the nonthermal gravitino production in 
supergravity, the most conservative lower bounds on the reheat- 
ing temperature, Tr > 3 — 4 MeV, comes from the insufficient 
thermalisation of the neutrino background, which changes the 
4 He abundance [2111 - 
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the end of inflation, and the onset of the field oscilla- 
tions respectively. We consider the two-stage dissipation 
mechanism described in the introduction and take 



C A T 



3 J.-2 



(8) 



where 2 C A = 0.64 /i 4 7V. 



A. Inflation and first reheating 

We are interested in the amount of expansion that fol- 
lows after the observable Universe exits the horizon, and 
so we require that warm inflation is supported for around 
50 e-foldings. In the warm inflation paradigm, the slow- 
roll equations that apply during inflation are 



V 



3H(1 + Q) 



and p r 



T0 2 



(9) 



These equations hold for as lon g as the so-called modified 
slow-roll parameters, given by [23| 



(10) 



(i + Q)' 



the kinetic density of the field, which in the strong dissi- 
pation regime remains subdominant during inflation, the 
Hubble parameter when the radiation density catches up 
with V((j}) is well approximated by 



H e ~ V2H fl' 2 



(15) 



Although the accelerated expansion finishes at this time, 
quasi-de Sitter inflation ends somewhat earlier. Never- 
theless, the amount of expansion that follows from that 
moment until the end of accelerated expansion is negligi- 
ble. Therefore, we approximate the number of e-foldings 
iV* by considering that the exponential expansion finishes 
at the time of potential-radiation equality. Also at this 
time, the source term Tip 2 in Eq. @ starts becoming 
subdominant, and the radiation density separates from 
its attractor solution and starts scaling as p r oc a~ 4 . On 
its part, and provided that Q e > 1, the scalar field still 
obeys its slow-roll equation and continues transforming 
its potential energy into radiation. Consequently, imme- 
diately after the time of equality the scalar density be- 
comes subdominant and the Universe becomes radiation 
dominated. We define the first reheating temperature Tr x 
as the temperature at the time of equality, namely 



T, 



Ri 



(16) 



V4> 



(i + Q) 



and 



1 V V 



(1 + Q) ZH 2 T 



(11) 



(12) 



are sufficiently small. In the above, and rj^ are the 
slow-roll parameters in the absence of a dissipative mech- 
anism 



m 2 p fV 



2 V V 



, ^ = m P 



v 



(13) 



Because the radiation density p r remains subdominant 
until the end of inflation, we may obtain the evolution of 
H during inflation simply by integrating d{H 2 ) ~ 5^2- ■ 
Using Eq. ©, the Hubble parameter when cosmological 
scales exit the horizon is 



H* ~ H f, 



1/2 



(14) 



where H is defined by Vq = 3H 2 m 2 p . In the strong 
dissipative regime of warm inflation, the accelerated ex- 
pansion lasts until the radiation density, constantly pro- 
duced owing to the dissipation mechanism, catches up 
with V ((/)). Somewhat before this time, the Hubble pa- 
rameter slightly deviates from H ~ Hof 1 / 2 . If we neglect 



2 This should be small in case of MSSM \22 



B. Onset of inflaton oscillations 

Because T depends on the temperature of the radia- 
tion, in the radiation dominated phase T evolves in a 
manner different than during inflation. As a result, the 
inflaton field loses energy at a different rate. Here we 
compute the potential density of the field when the field 
ceases to follow its attractor equation. This time, owing 
to the dissipation term, is determined by the condition 



T ~ |U 'f/ 2 



(17) 



Integrating the slow-roll equation </> ~ |V'|/T from the 
end of inflation until the onset of the fast-roll motion we 
obtain the integral equation that determines <p 



V'(p 2 



dt 



= It, 



(18) 



where we used T = C^T 3 /^ 2 . Using that in this stage 
the Universe is radiation dominated, i.e. H ~ H e (a e /a) 2 
and T ~ Tfj 1 (a e /a), the integral It may be approximated 
by 



da 



1 



{aH)C^T 3 



1 



(19) 

The ratio a a /a e can be estimated as follows. Using the 
definition of T and that T a a -1 it is straightforward to 
obtain the relation 

To-Tef^Vg^, (20) 
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where T e can obtained from the attractor equations and 
T Q is determined by Eq. (|17p . The redshift factor (a /a e ) 
is then approximated by 

5 C 20/21 (A)10/7(1 _ ^ ) 10/7 if 10/21 m 25/21 

' 2a 5/7 i/"ol 5/6 /e 5/14 (i-5 o ) i0 / 3 ^ /3 ' 

(21) 

where ' = -4. Substituting this expression into Eq. (fT9|) 
. solve Eq. (TT8")) numerically and find 5 a . 



Oo 



we can 



C. Second reheating temperature 

In order to simplify our analytical treatment we assume 
in this section that at a = a a the field starts performing 
fast oscillations about its vev. In such case its average 
density (p$) can be readily computed (24|. Although this 
assumption naturally results in an error estimating the 
scalar density of the field, in the Appendix we explain 
how this can be corrected, thus obtaining an accurate 
estimate of the average scalar density at late times. As- 
suming then a sudden transition to the stage of fast os- 
cillations, the average density of the field is determined 
by 



(p <t> )^-3H(l + Q)(p 4> ). 



(22) 



Although at the beginning of the oscillations it may be 
Q > 1 , we keep the term 3H in the above as this becomes 
the dominant one when the system moves into the weak 
dissipation regime. Using that H ~ H Q {—) and that 

((f)) = 4> v , and hence T ~ T Q (^f) 3 , the average density 
of the field can be readily obtained 

( P0 )^(^) o (^) 3 exp{-3Q o (l-^)} , (23) 

where Q can be expressed in terms of 5 e and <5* after 
using Eq. (|20l) . From the above expression it follows that 
(ptj)) oc a -3 a few Hubble times after the onset of oscilla- 
tions. Thus, provided that Q Q is not large, the inflaton 
may come to dominate the Universe at a later time well 
before nucleosynthesis. We denote by i eq the time when 
the average scalar density (p^) catches up with p r . The 
subsequent matter-dominated phase must be sufficiently 
long so that the entropy produced when the inflaton field 
decays completely at t = tdec is enough to dilute the over- 
production of gravitinos. With the complete decay of the 
inflaton at t = idee, which we consider a free parameter, 
the Hot Big Bang evolution is recovered. The second 
reheating temperature Tr 2 is then 



Ado 



3/4 



(24) 



which depends on S e and on the model parameters. Obvi- 
ously, to estimate this temperature we first need to have 
an estimate for the quantities (5* and 6 e . 



III. OBSERVATIONAL CONSTRAINTS 

The purpose of this section is to determine (5* and 8 e , 
employed to compute Tr 2 , in terms of the model param- 
eters. The raison d'etre of this is to link the reheating 
temperatures Tr 2 , computed according to Eq. ((24"]) . to 
observable quantities like the spectral index or its run- 
ning, in turn determined by the model parameters. By 
doing this we manage to identify the interval of tempera- 
tures Tr 2 corresponding to certain range of values of the 
spectral index (see Fig. [5]). 

The predicted amplitude of the curvature perturbation 
spectrum in the strong dissipativc regime of warm infla- 
tion is given by 



where 



2ttV V H 



(25) 



(26) 



has been known for sometime [loj . whereas the correc- 
tion was computed only recently [25[ . The correction 
term has a range of behavior, depending on the details of 
the dynamics of the radiation energy density component. 
We have examined this range, but will not explore these 
details here. The purpose of this paper is to highlight 
a general mechanism for treating gravitino production 
within warm inflation, and for that we will focus on one 
particular form of the correction [25j, |26| for the dissipa- 
tion coefficient in Eq. 



(c) 
K 



1 + a Q 2 



<2l( 



(27) 



s/T+Q 

where a ~ 0.662 and a\ ~ 3.26 x 10 -4 . 

Using Eqs. (j6]), ((!]) and (fT4|) . and writing Q in terms 
of C,/,, Hq and <f> the above may be recast as 



(28) 



which can be solved numerically to determine </>*. Here 
/3j denotes the rest of the model parameters which the 
amplitude of the spectrum may depend on. 

In the strong dissipative regime warm inflation gives 
way to a radiation dominated Universe. Neglecting the 
kinetic density of the field, the end of the accelerated 
expansion occurs when the radiation density catches up 
with the scalar potential, i.e. when 



V e = (p r )e 



(29) 



As mentioned earlier, the amount of quasi-de Sitter infla- 
tion is well approximated by considering that this finishes 
at 4> = 4> e . To determine <f> e we rewrite the above condi- 
tion expressing V e in terms of S e and using the attractor 
equation for (p r )- The field value 4> e is then determined 
by 

(A) 8 / 9 (l-^e) 8 / 9 ^ 2 11/9 ^ /9 (Ho\ 2/9 
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By simple inspection, it is clear that solving for 0* and 
4> e from Eqs. ([25)1 and ([50)1 is, in general, not amenable 
to analytical treatment. Owing to this, we proceed al- 
ternatively by solving Hq and in terms of (5* and <5 e , 
i.e. 



#o = H (6*,6 e ) and 0^ = 0^(8^,8^. 



(31) 



The amount of exponential expansion that follows from 
horizon crossing is given by N ~ J -pPj dcj). Although this 
approximation is valid in the strong dissipative regime 
(i.e. Q 3> 1), by solving numerically the inflationary 
stage we find that this formula still constitutes a good 
approximation to the number of e-foldings for > 2. 
Writing H ~ Hq J 1 / 2 (which holds until the end of infla- 
tion) and using Eq. (f5CT|) . the amount of expansion after 
the observable Universe exits the horizon is given by 



(A) 8 / 7 (l - s e f 7 
2 i6/7 f e n 



5, 



f/ 7 dS 



Se (1 - 5)8/7(/')V7 



(32) 



Keeping fixed the number of e-foldings N* it is possible 
to solve for <5* in terms of 6 e by integrating numerically 
the above equation. Our strategy now consists of em- 
ploying this parametrisation <5* = <5*(<5 e ) and Eq. (|3ip to 
express the observable quantities in terms of S e . 

In the absence of the corrected part of the pcrturba- 

(c) 

tion spectrum, i.e. V-j^ = 1, the spectral index of the 
perturbation spectrum, defined as n s — 1 = cU J\ v 7 l , is 
given by [23j 



ra. - 1 



where 



(-(2 - 5A)e - SA^ + (2 + 4A)ct ) , 

(33) 



'Vcj) 



and A 



Q 



1 + 7Q 



(34) 



Using the slow-roll equations, the ratio Q at the time of 
horizon crossing is given by may be rewritten as follows 



Q * 



2 6/7 3 l/7( 1 _ ^)8/7 j5/7 \ \ mp 



(35) 

Upon including the correction Vj^ , the spectral index 
of the corrected spectrum is given by 



n s =n s + 



dlnV- 



(c) 
7? 



din k 



(36) 



where 



-H- A(lOe - 6r U + 8a ) . (37) 



din k 



dQ 



Using the combined WMAP+BAO+SN data and for 
negligible tensor perturbations, current observations con- 
strain the spectral index to n = 0.963 ± 0.014 at the 1-cr 
level [27 1 . 



Also interesting for observational purposes is the run- 
ning of the spectral index, defined as 



n a = n. +n 



'(c) 



(38) 



where nj can be written in terms of slow- roll parameters 
0] and 



Ac) 



dlnV- 



d 2 \nV^ ] d 2 lnP, 



(c) 



dlnfc 2 



dQ- 



2- (lOe - 6^ + 8ct ) 2 A 2 



(c) 
■R 



dQ 



(lO6 -6?? + 8cr ) S 



(1 + 7Q) 2 



+ 10e^, - 6770 + 8cr0 



-4. 

(39) 

In the absence of tensor perturbations and using the 
combined WMAP+BAO+SN data, the running is con- 
strained to the interval n' s = —0.034 + 0.026 at the 1-a 
level (13. 



IV. CONSTRAINTS ON GRAVITINO 
PRODUCTION 

Here we discuss production processes of gravitino and 
observational constraints on its abundance. Fist of all, 
the gravitino is produced through the thermal scattering 
among the standard particles such as quarks and gluons. 
Then the yield value, Y 3 / 2 = n 3 / 2 /s of the thermally- 
produced gravitino is approximately given by [2814321 ] 



Y ^ * 2 x 10 " 4 (lo^v) 



3l7l l/2, 



(40) 



with rrig being the gluino mass. Note that this yield value 
is proportional to the reheating temperature. For mas- 
sive unstable gravitinos with m 3 / 2 3> 10 2 GeV ~ £)(m§), 
the second term in the second bracket is negligible. The 
energetic products by the decay of gravitino are so dan- 
gerous that the light element abundances are modified 
by their annihilations and productions through scatter- 
ing processes. When we consider the case of m 3 / 2 ~ 0(1) 
TeV, to agree with the observational light element abun- 
dances we get upper bounds on the gravitino abundance 
to be Y 3/2 < 10- 16 and Y 3/2 < 10" 14 for a branching 
ratio deca ying into hadrons, Bh = 1 and Bh = 10 -3 , re- 
spectively |29l l33j. The latter value of Bh is a reasonable 
lower limit on the hadronic branching ratio [29| . Then 
we obtain the upper bound on the reheating tempera- 
ture in turn through Eq. (gDJ) to be T R < O(10 6 ) GeV 
(T < O(10 8 ) GeV) with a hadronic branching ratio 
B h = 1 [B h = 10-3) [U, HI (see Fig. 44 and Fig. 45 
ofRcf. [2g|). 

On the other hand, for much larger masses, m 3 / 2 ^ 70 
TeV with B h = 1 (or m 3/2 > 20TcV with B h = 10~ 3 ), 
we have another type of upper bound on the abundance 
of gravitino which produces the Lightest SUSY Particle 
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(LSP) dark matter as a decay product. This is estimated 
to be 



Y 3/2 = Y LSP < 4 x 10- 



_12 {QlspIi 2 \ ( m LS p \ — 1 



)( 



0.12 J V10 2 GcV 



.(41) 



By adopting the upper bound on the density param- 
eter of LSP, n L spli 2 < 0.12 with h the reduced 
Hubble parameter, and a typical scale of the LSP mass 
(tolsp ~ 10 2 GeV), this gives us the upper bound on the 
reheating temperature T R < 2 x 10 10 GeV. We see that 
this is much milder than that of BBN for the thermally- 
produced gravitinos with m 3 / 2 ~ 1 TcV. 

Other attractive process to produce gravitino would 
be the nonthermal production by the inflation decay 
through 4> — > 2^ [3J-[44| where ip^ means the grav- 
itino. Since we expect a relatively low reheating tem- 
perature for the second reheating, this process can im- 
portantly affect the result because the abundance of 
the nonthermally-produced gravitino is inversely propor- 
tional to the reheating temperature, Y3/2 ~ 2B^/ 2 Y^, ~ 
3 .A 1 pT '0^2V> M / '2m^Tfl where r^-^i/v is the differential de- 
cay width to a pair of gravitino, -B3/2 = T^-y^ /T^ is 
the branching ratio into a pair of gravitinos, and is the 
total decay rate of the inflaton filed (~ T' R /M p ). Adopt- 



ing a result of Ref. [44|, the yield of the nonthermally- 



produced gravitino is estimated to be 



3/2 



- 7 x 10~ 15 
( T R V 

Vio 8 Gevy 



(- 



-V 



10 18 GeV/ V10 8 GeV/ 



(42) 



in case of < ^jm 3 / 2 m p where (<fi) and are the vev 
and the mass of the inflaton at the decay epoch. There- 
fore the yield value explicitly depends on the model pa- 
rameters such as (</>) and m^, which means that we need 
to analyse the effect for every models. We will also dis- 
cuss this type of the constraint on some models simulta- 
neously in the following sections. 



V. SOME EXAMPLES 

In this section we illustrate the procedure developed 
in the previous section in order to compute the second 
reheating temperature Tp 2 in some simple models. 



A. Tree level potential 

We examine a simple extension of the hilltop model 
studied in Ref. [2(| where we consider an additional quar- 
tic term 



V = V - 1 -\m 2 \tf + ^\ 



(43) 



that stabilises the potential and provides a minimum to 
support the field oscillations after inflation. The param- 
eter A is tuned so that the scalar potential vanishes at 
the minimum. With this choice, the scalar potential is 
written as in Eq. (|6]) with 



/ = S 2 (2 - S) 2 



(44) 



and defining 1 770 j = |?ri 2 |/3i?| 2 the expectation value of 
the field is 4> v = 2\r) \~ 1 / 2 mp. 

In order to display the typical range of parameters al- 
lowed by observations (see Fig. [T]) we fix the amount of 
inflation after horizon exit to the typical value iV* ~ 50. 
Using the parametrisation 5* = <5*(<5 e ) obtained from 
Eq. (j32|) , the spectral index and its running are deter- 
mined by S e and 770 . We use Eq. ([30]) to plot our results 

in Fig. Q] in the plane ?7o-log 10 C^,(Vq^ 4 /mp). In the 1-er 
window of n s displayed in Fig. [T] we find 



10 1U GeV < H < 10 12 GeV 



(45) 



which corresponds to the range 10 7 < < 10 s and 
10 8 GeV < to < 10 10 GeV. Using Eq. ([16]) to compute 
the first reheating temperature Tp x , we find that this 
ranges between 



10 13 GeV<T ff , < 10 14 GeV. 



(46) 



This is then far in excess of the current bounds for ther- 
mal gravitino production [2^, [33|. Also, in the entire 
space depicted in Fig.Q]the running of the spectral index 
ranges between — 10 -3 < n' s < — 10 -4 . 



Before estimating Tp 2 by using Eq. (|24|) . we need to 
obtain an accurate estimate of the average scalar density 
(ptf>) at late times. As mentioned before, the estimated 
Tp 2 in Eq. (|24|) is not accurate because of two reasons. 
The first one owes to have assumed that the field un- 
dergoes a sudden transition to the phase of fast oscilla- 
tions, whereas the second one arises because the condi- 
tion that determines the onset of the fast-roll motion it- 
self involves an order of magnitude estimate [cf. Eq. ([T7|) ] . 
As a result, the redshift factor (a /a e ) in Eq. (f2Tj) can- 
not be computed accurately. Using Eq. ([2U|) and writing 
H Q ~ H e (a e /a ) 2 we obtain that Q Q is given by 



Qo — Q< 



(l-Se) 2 fa, 
(1~S ) 2 U 



(47) 



Consequently, the ratio Q can be determined only up to 
a factor of order unity, which can change the reheating 
temperature Tp 2 substantially owing to its exponential 
dependence on Q a [cf. Eq. ([23])]. Note also that although 
5 , obtained by solving numerically Eq. (|18j). cannot be 
accurately computed either, this docs not imply a signif- 
icant error in Q Q because S Q <C 1 in any case. 
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FIG. 1: Region of the parameter space allowed by observa- 
tions. Heavy-shaded areas are excluded either because the 
inflaton field decays completely after inflation (Q e > fO) or 
because the system is in the weak dissipation regime (Q < 1). 
The light-shaded area identifies the region where the inflaton 
field does not decay completely after inflation (I < Q e < 10) 
and the spectral index is outside its 95% CL interval. In the 
unshaded region the inflaton avoids its complete decay after 
inflation and the spectral index is within 95% CL. Dashed 
lines (red) correspond to constant values of the ratio Q„, 
thin lines (blue) correspond to constant values of n s and dot- 
dashed lines (green) correspond to constant Hq. 




0.90 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 



n s 

FIG. 2: Reheating temperature Tr 2 vs spectral index n 3 for 
several values of | t?o | - We set set m 3 / 2 ~ TeV and Bh — 1CP 3 
(see Fig. [3]). Only the part of the curves compatible with 
strong dissipation is shown. The horizontal shaded band cor- 
responds to the range where the problem of thermal overpro- 
duction of gravitinos is less severe pffll . |33| ] The vertical band 
corresponds to the preferred range of the spectral index at 
l-a. 



-4 -2 2 4 




-4 -2 2 4 
Log 10 (m 3/2 /GeV) 



-4 -2 2 4 




Log 10 (m 3/2 /GeV) 



FIG. 3: Plot of A CI1 t = T C q/Tn 2 for a hadronic branching 
ratio, (a) Bh = 1 and (b) Bh = 1CP 3 , which is necessary 
for sufficient dilution of the gravitino abundance after the 
second reheating to agree with the observational constraints. 
From the top to the bottom, we plot the cases of Tr 1 = 
10 14 , 10 12 , 10 10 and I0 8 GeV, respectively. 

In the Appendix we show that in order to obtain an 
accurate estimate for Tr 2 it suffices to perform the trans- 
formation 

Q -> Q = aQ , (48) 
where a ~ 1.7. Applying such a transformation, we use 
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the late-time corrected densities in Eqs. (|6Tj) and ([62]) to 
estimate Tr 2 from Eq. (|5H 

The decay of the infiaton field must produce enough en- 
tropy as to dilute the excess of thermal gravitinos pro- 
duced at temperature Tr x . Given a certain reheating 
temperature T Rl and for a certain gravitino mass rfi 3 / 2 , 
this is achieved provided the second reheating temper- 
ature is not larger than the upper bound discussed in 
Sec. [GO 

T R2 = A ent( m 3/2, T Rl )T cq . (50) 

Here A en t is the entropy dilution factor and T oq is the 
temperature of the radiation bath at the time of matter- 
radiation equality, which is the same order of Tr 1 . In 
Fig. [5] we present our results for Tr 2 vs n s for several 
values of 1 770 1 and for the typical values JV» ~ 50. To ob- 
tain the curves in Fig. [5] we set m 3 / 2 ~ O(TeV) and 
B h = 10~ 3 , hence A ont - 10 6 (T Kl /10 14 GeV). Further- 
more, we arrange the complete decay of the infiaton at 
the earliest time compatible with the dilution of the grav- 
itino overproduction, i.e. we set (a eq /adcc) 3 ^ 4 = A^ 
Eqs. JH]) and (HO])]. If the decay of the infiaton is fur- 
ther delayed then Tr 2 decreases accordingly. Note also 
that an excessively low reheating temperature may pose a 
problem when the non-thermal production of gravitinos 
is considered, therefore it is desirable that the infiaton 
field decays "soon" enough to obtain a sufficiently large 
reheating temperature. 

The shaded areas displayed in Fig. [5] enclose the 1- 
a window of n s and the range of temperatures where 
the thermal gravitino overproduction are typically less 
problematic, i.e. Tr 2 ~ 10 6 GeV for Bh = 1 and 
TR 2 ~10 8 GeV for B h = 10~ 3 [H, HI- 0ur P lots 
make clear that, provided 1 770 1 lies in the range 0.02 < 
\Vo\ ;$ 0.05, the infiaton, while giving rise to a ther- 
mal perturbation spectrum in agreement with current 
observations, manages to drive a late matter-dominated 
epoch and give rise to a radiation bath with temperature 
Tr 2 < 10 6 — 10 s GeV after its complete decay. This is an 
interesting result, and hence it is worth mentioning that 
the "survival" of the infiaton field owes to the moderate 
growth of Q during inflation. For the hilltop potential 
in Eq. (|43p . in the slow-roll regime, the dissipative ratio 
evolves as Q: 

1n7 TTtq y ~ \mp~) ) ' (51) 

and thus Q starts decreasing and increases only towards 
the end of inflation when cf>/mp ~ 1 770 1 1 / 2 . In particu- 
lar, we find 1 < Q e /Q* < 3.5 within the allowed region 
displayed in Fig. [TJ Ultimately, the fact that T hardly 
grows during inflation is due to the moderate steepness 



of the potential. Therefore, to avoid the complete de- 
cay of the infiaton it is only necessary to tune the model 
parameters so that the last stage of inflation takes place 
with the system not far away from the weak dissipation 
limit, and therefore for not too large values of [ 770 1 - We 
note however that for the range 0.05 < 1 770 1 < 0.10 where 
the gravitino overproduction is avoided, the vacuum ex- 
pectation value of the field <f> v = 2\rjo\~ 1 ^ 2 mp is slightly 
above the Planck scale. 

In Fig. [3] A on t is plotted vs m 3 / 2 for several values of 
T Rl for (a) B h = 1 and (b) B h = lO" 3 . For m 3/2 > 
10 2 GeV, we adopted the constraints from the unstable 
gravitinos [2f|. On the other hand, for m 3 / 2 < 10 2 GeV, 
we adopted the upper bound on the density of the stable 
gravitino not to exceed the LSP density shown in Eq. (|4"Tj) 
by replacing mLSP to m 3 / 2 , and using Eq. (|40|) . Then the 
second term in the second bracket of Eq. flU)) dominates. 
The left-side tail with respect to the peak structure in 
Fig. E] comes from the saturation of the number density 
due to the thermalisation of gravitino, V 3 / 2 < 1/g* . 
In the current cases where T Rl < 10 14 GeV and 777.3/3 = 
O(TeV), the required dilution factors are of the order of 
Ac„t > O(10 8 ) and > O(10 6 ) for B h = 1 and 10~ 3 , which 
implies the constraint Tr 2 < O(10 6 ) GeV and Tr 2 < 
O(10 8 ) GeV, respectively. 

From (j43]) and m ~ ^/\r]o\Ho, we find that the infia- 
ton masses would be in the range of 10 8 GeV < m<p < 
10 10 GeV because ra^ ~ m at the oscillating epoch. By 
considering the constraint from the nonthermal produc- 
tion of the gravitino given in Eq. (|42j) with (0) ~ m p 
and ni (p - 10 s GeV, we get Tr 2 > 10 8 GeV (Tr 2 > 10 10 
GeV) for B h = 10 -3 (B h = 1). It is attractive that the 
gravitino abundances produced both thermally and non- 
thermally agree with any observational constraints if we 
adopt the hadronic branching ratio to be Bh = 10~ 3 . 

The strong dissipation limit of warm inflation is known 
to result in the generation of Non- Gaussian effects (45I — 
\at\ , and several models of warm inflation have been con- 
structed with such effects @, 5H . In such a limit it 
was shown in Ref. [4H ] that entropy fluctuations during 
warm inflation play an important role in generating non- 
Gaussianity, with the prediction 

_ 15h ( 1+ ^< /NL <| ta ( 1+ !)4 

(52) 

To estimate the magnitude of /nl for this model within 
the allowed region shown in Fig. [2] it is enough to 
compare the results plotted there (left-hand panel) with 
Fig. [2 We then see that in the allowed region (J* is 
at most of order 10, which corresponds to |/nl| < 10. 
This is well within the observed range of the /nl = 
32 ± 21 (68%C.L.) [13. Note that the PLANCK satel- 
lite [isl , which was launched recently, will be sensitive to 
non-Gaussianity at the level |A/ NL | = 0(5) |49H50j. 



B. Supergravity inspired model 

We consider now a model typical of supergravity the- 
ories [U as an example. In this kind of models, 
the inflaton supcrficld $ is assumed to have an R charge 
2/(n + 1) allowing the superpotential 



Wo 
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1 



(53) 



with n positive and g a coupling constant. The con- 
tinuous U(1)r symmetry is assumed to be dynamically 
broken to a discrete Z 2n R at a scale v -C nip generating 
the superpotential 



ll- 



w 2 $ - 



1 



n+1 



(54) 



The inflaton field c/>(x)/y/2 is then identified as the real 
part of the scalar component of the supcrficld <£>. Taking 
the i?-invariant Kahler potential, K = |$| 2 + fc|$| 4 /4 + 
• ■ ■ , the scalar potential is given by 



v{4>) 



g 2 4> 



2)i 



2 n/2- 



(55) 

which in the absence of dissipation was shown to be flat 
enough to support inflation and generate the appropriate 
perturbation spectrum. At the vacuum, one has 



\ m p9, 



mp . 



(56) 



and the scalar potential is negative due to the contribu- 
tion proportional to k from the non-minimal Kahler po- 
tential. In the original model [5l| that was canceled by 
a positive SUSY breaking effect contributing an amount 
Ag USY , which then fixes the gravitino mass. Thus, the 
main role of the non-minimal Kahler k term was to fix 
the scale for the gravitino mass. On the other hand, 
by considering k ~ g and v <C mp, one has in general 
v -C <f> v , and the quadratic term in Eq. (|55)1 can be ne- 
glected whenever <f> ~ (f> v . This is what we expect during 
the last 50-60 e-folds of warm inflation. As discussed 
earlier, warm inflation finishes giving way to a radiation 
dominated Universe when <f> < <p v , while the field is still 
in slow-roll. Owing to the slow-roll motion, it is natu- 
ral to expect that v -C (f> < <j> v . Consequently, the scalar 
potential during the last stage of inflation can be approx- 
imated by neglecting the quadratic contribution due to 
the coupling k. 

In the following we will just set k — to discuss this 
kind of scalar potentials in the context of warm inflation. 
This also means that we do not consider any particu- 
lar susy breaking mechanism, neither we link the grav- 
itino mass to any particular vacuum scale. The potential 

1 and 



Eq. ([55)1 can be written as in Eq. © with Vq 



(57) 



N,=50 




0.05 0.1 0.15 

3g(m p /v) 2 /2 1 ' 2 

FIG. 4: Parameter space available to the strong dissipative 
regime of warm inflation for the sugra potential Eq. (|55[) . with 
n = 3 and k — 0. In the left shaded area the system is in the 
weak dissipative regime (Q e < 1), while the right shaded area 
is excluded because the inflation field decays completely after 
inflation (Q e > 10). The dashed lines (red) correspond to 
constant values of Q*, while the thin (blue) lines are those of 
constant n 3 . The short-dashed line gives that of constant field 
end value <f> e - Below that line the field is always subplanckian, 
and it is mp < <j> < Wrap above. The dot-dashed (green) 
lines are those of constant vacuum energy V 1//4 = v. 



and we can apply the procedure developed in the last 
section to find the range of model parameters consistent 
with observations. For a general power n, the dissipative 
ratio evolves during inflation as: 



dlnQ 

dN R 



\m\/(n-i) 
1 + 7Q 



14- 6n-5 



n — I \ nip 



(58) 

When n = 2 we recover the potential studied in the pre- 
vious subsection, with \t)q\ = gv 2 /m P , and a moderate 
decrease/increase of Q during inflation. On the other 
hand, for steeper potentials with n > 2, Q always in- 
creases. Even if the observable universe exits the horizon 
when the system is still in the weak dissipative regime 
(i.e. Q* ^ 1): inflation typically finishes into the strong 
dissipative regime. We are interested here on finding 
the model parameters which lead to a moderate value 
Q e - O(10) when n > 2. 



Fig. U depicts the range of parameters allowed by ob- 
servations for n = 3 and ~ 50. In this case the in- 
crease of the dissipative ratio Q is still moderate and we 
can find values of parameters which leads to the strong 
dissipative regime but with Q e < 10, and for which the 
prediction for the spectral index is within the observa- 
tional range. For larger powers n > 4, owing to the 
steepness of the potential, a large ratio Q is needed to 
drive slow-roll inflation. As a result, a substantial part 



10 



of the region where the spectral index agrees with ob- 
servations would be excluded by the bound /nl < 74 
[I2}, implying < 1.2 x 10 4 . And in regards to the 
gravitino overproduction, the substantial increase of the 
ratio Q during the inflationary stage reduces consider- 
ably the parameter space for which Q e ~ O(10). For 
most of the parameter space, inflation finishes well into 
the strong dissipative regime with Q e 1. But even if 
we tune the parameters to avoid too large a value of Q e , 
the model with n > 4 predicts a blue spectrum and no 
tensors. This is in conflict with observations, which in the 
absence of primordial tensor perturbations favour a red 
tilted spectrum. Nevertheless, we remark that the blue 
spectrum stems from the large correction to the per- 
turbation spectrum. If one disregards such a correction, 

(c) 

i.e. = 1, the model gives rise to a red tilted spec- 

trum with 0.956 < n s < 0.974, in agreement with current 
observations. Still, the dissipative mechanism with too 
step potentials does not help with the gravitino problem. 

VI. CONCLUSIONS 

In this paper we have shown that warm inflation mod- 
els can lead to a new mechanism for controlling gravitino 
overproduction. The residual oscillating energy of the 
inflaton field after inflation eventually dominates the en- 
ergy density of the universe. Then the late-time entropy 
production by its decay can really dilute the gravitinos 
which are thermally-produced in the radiation dominated 
epoch just after warm inflation ends. Even if we consider 
the nonthermal production of gravitino, this scenario is 
consistent with the observational constraints. We have 
demonstrated that this mechanism is applicable to a large 
class of models, when the dominant term during inflation 
goes like 4> n with n — 2 , 3, with a mild decrease/increase 
of the dissipative parameter during inflation such that 
still Q e < 10 by the end of inflaton. For steeper po- 
tentials, the increase of the dissipative parameter would 
lead to Q e 3> 10 for most of the parameter space, and 
the complete decay of the inflaton by the end. 

We have also discussed the possibility to detect the 
non-Gaussianity of the order of /nl ~ 10 which origi- 
nates from the strong dissipation in the current models 
of warm inflation. The PLANCK satellite will be able 
to detect this signature by which we can distinguish the 
current model from the normal cold inflation models. 
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VII. APPENDIX 

In this section we obtain an accurate estimate for Tr 2 . 
In order to do so, we need a reference time to compare 
the predicted scalar density after inflation, Eq. ([2"3")l , with 
the value obtained from the numerical solution of the 
system, Eqs. |T]) and @. An appropriate "checkpoint" 
is provided by the time when the system moves into the 
weak dissipation regime: T w = 3H W (or Q w = 1). This 
checkpoint is appropriate because when T = 3H the field 
is already performing fast oscillations, hence the density 
is well approximated by the average in Eq. (|2"3"|) . Also, 
for the cases of interest (when the inflaton does not decay 
completely right after inflation), the system is not far 
away from the weak dissipation regime. This then allows 
us to obtain numerically the scalar density at this time 
quite easily, as T = 3H not too late after the field starts 
oscillating. 

Using Eq. (|2"3"|) , the predicted average scalar density 
when the system reaches the weak dissipation regime is 

( P4> ) w ~ (pt) Qo 3 exp {-3 Q (1 - Q- 1 ) } , (59) 

where a Q /a w — Q^ 1 since Q falls as a -1 after inflation. 
Because the ratio Q a is determined up to a factor of order 
1, it is possible to match the above prediction to the 
numerical solution at the time a = a w by the performing 
the substitution 

Q ^Qo = aQ , (60) 

and then finding an appropriate value for a. Within the 
allowed region displayed in Fig.[T]and fixing (p<p) = 2U Q , 
we find that the average scalar density is matched to 
the numerical solution by taking a ~ 1.70. Hence, at 
late times (a 3> a w ) the average scalar density and the 
radiation density are well approximated by [cf. Eq. (|23l) ] 

( P<p ) ~ (pMQo)- 3 (^) 3 cxp{-3Q } . (61) 

and 

Pr ^2( Pr ) e (Q )- 4 ^ 4 (^) 4 , (62) 

with the redshift factor (a e /a ) as given by Eq. (|2ip . 
The factor of 2 in the last equation is necessary so that 
p r matches its numerical solution at a = a w . Although 
put by hand, the introduction of such a factor is justi- 
fied because after potential-radiation equality most of the 
scalar density, i.e. an amount V e = (p r )e, is transferred 
to the radiation bath in a Hubble time or so. 
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